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ABSTRACT 


Equilibria of complexation of Co(ID, Ni(ID) and Cu(ID ions with L-aspartic acid and ethylenediamine were investigated pH metrically in 0.0-60.0% v/v propylene 
glycol-water mixtures maintaining an ionic strength of 0.16 mol L-1 at 303.0+0.1 K. Stability constants of ternary complexes were calculated and various models were 
refined with MINIQUAD75. The best fit chemical models were selected based on statistical parameters and residual analysis. The species detected were MLXH, MLX 
and ML2X for Co(II), Ni(II) and Cu(II). The appropriateness of experimental conditions was verified by introducing errors intentionally in the concentrations of 
ingredients. The stabilities of the complexes followed the Irving- Williams order, i.e., Co(II) < Ni(II) << Cu(II). 
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1. INTRODUCTION: 

Several studies on ternary stability constants of a-amino acids have been 
reported in different media.'’Metal ions exist in non-exchangeable form, loosely 
bound to some biological ligands and in equilibrium with a variety of bioligands 
in different biofluids like blood serum, intestinal fluid, cerebrospinal fluid, gas- 
tric juice etc. Thus, simultaneous equilibria involving a variety of metal ions and 
ligands is possible in bio fluids. Cobalt, nickel and copper have several biological 
functions.”’ Therefore, investigation of ternary complexes containing a metal ion 
and two different ligands drew the attention of several researchers. Mixed ligand 
complexes can be considered as models for apoenzyme-metal ion-substrate com- 
plexes. In peptides, proteins and amino acids, amide group is one of the impor- 
tant binding sites for the coordination of the metal ions. Acidity and basicity ofa 
molecule is governed by its structure and solvent effects.’ Propylene glycol (PG) 
is chosen as a polar organic solvent to mimic the permittivity of the biological flu- 
ids. L- Aspartic acid (Asp) and ethylenediamine (en) are chosen as model com- 
pounds to proteins and substrates. 


Asp, a non-essential amino acid, plays an important role in maintaining the solu- 

bility and ionic character of proteins.° en is used as an important monodentate, 

bidentate or a bridging ligand.’ It is useful in manufacturing accelerator or curing 

agent in epoxy industry. It is involved in the synthesis of B-enaminoesters’ and 

Schiff bases.’ The protonation constants of ethylenediamine were reported ear- 

lier by theoretical calculations.'*"' Protonation constants and binary stability con- 
stants of Asp" and en” with Co(II), Ni(II) and Cu(II) in PG-water mixtures were 

reported earlier. Hence, chemical speciation of their ternary complexes is 

reported in this communication. 


2. EXPERIMENTAL: 

Aqueous solutions (0.1 mol L") of Co(II), Ni(II) and Cu(II) chlorides (GR Grade, 
E-Merck, Germany) were prepared by dissolving them in triple distilled water. 
0.05 mol L" aqueous solutions of L-aspartic acid (GR Grade E-Merck, Germany) 
and ethylenediamine (AR, Qualigen, India) were also prepared. To increase the 
solubility of the ligands and metal salts, 0.05 mol L" hydrochloric acid was main- 
tained in the solutions. Propylene glycol (Finar, India) was used as received. The 
strength of acid in metal ion solutions was determined using the Gran plot 
method.'*"’Errors in the concentrations of the ligands, metal ions and alkali were 
subjected to analysis of variance (ANOVA)..” 


Titrations were carried out in the medium containing varying concentrations of 
PG maintaining an ionic strength of 0.16 mol L" with sodium chloride at 
303.0+0.1 K, with ELICO (Model LI-120) pH meter of 0.01 readability in con- 
junction with a glass and calomel electrode. The pH meter was calibrated with 
0.05 mol L" potassium hydrogen phthalate in acidic region and 0.01 mol L borax 
solution in basic region. The glass electrode was equilibrated in a well stirred PG- 
water mixtures containing inert electrolyte. The effect of variations in asymme- 
try potential, liquid junction potential, activity coefficient, sodium ion error and 
dissolved carbon dioxide on the response of glass electrode were accounted for in 
the form of correction factor (log F) which was computed from the experimental 
and simulated acid-base titration data calculated by SCPHD program." A correc- 
tion was applied to the pH meter dial readings to accounts for the solvent effect 
on pH. 


Titration of strong acid with alkali was carried out at regular intervals to check 
whether complete equilibration was achieved. The calomel electrode was 


refilled with PG-water mixtures of equivalent composition as that of the titrand. 
In each of the titrations, the titrand consisted of 1 mmol of hydrochloric acid in a 
total volume of 50 mL. Titrations were carried out in the presence of different rela- 
tive concentrations of the metal (M) to Asp (L) to en (X) (M:L:X = 1:2.5:2.5, 
1:2.5:5.0, 1:5.0:2.5) with 0.4 mol L’ NaOH (Table 1). The best-fit chemical 
model for each system investigated was arrived at using MINIQUAD75”, which 
exploits the advantage of constrained least squares method in the initial refine- 
ment and reliable convergence of undamped, unconstrained Marquardt algo- 
rithm. The protonation constants and the stability constants of the binary metal 
complexes of the ligands were fixed in refining ternary complexes. 


3. RESULTS AND DISCUSSION: 

3.1. Modeling Strategy: 

A preliminary investigation of alkalimetric titrations of mixtures containing dif- 
ferent mole ratios of Asp and en in the presence of hydrochloric acid and inert 
electrolyte indicates that no condensed species were formed. The existence of the 
ternary complexes was determined by performing an exhaustive modeling study 
and the results of a typical system are given in Table 2. The models were evalu- 
ated assuming the simultaneous existence of different combinations of species. 
Models containing various number and combinations of species were generated 
using an expert system package CEES and these models were refined using 
MINIQUAD7S. As the number of species increased, the models gave better sta- 
tistics denoting the best fit. The best fit models were chosen based on the statisti- 
cal parameters like y’, R-factor, skewness and kurtosis given in Table 3. The ter- 
nary complex species detected are MLX, ML,X and MLXH for Co(II), Ni(ID) 
and Cu(II). The present values are compared with those reported earlier studies 
(Table 4). 


Avery low standard deviation (SD) in the overall stability constants (log B) indi- 
cates the precision of the parameters. The small values of U,,,, (sum of squares of 
deviations in the concentrations of the metal, the ligands and the hydrogen ion at 
all experimental points corrected for degrees of freedom) indicate that the mod- 
els represent the experimental data. Small values of mean, standard deviation 
and mean deviation for the systems corroborate that the residuals are around a 
zero mean with little dispersion. For an ideal normal distribution, the values of 
kurtosis and skewness should be three and zero, respectively. The kurtosis values 
in the present study indicate that most of the residuals are very nearer to 
leptokurtic and a few form mesokurtic patterns. The values of skewness recorded 
in Table 3 are between -0.46 and 3.46. These data evince that the residuals form a 
part of normal distribution hence, least-squares method can be applied to the 
present data. The sufficiency of the model is further evident from the low crystal- 
lographic R-values recorded. 


3.2. Solvent effect on stability of ternary complexes: 

PG isa protophilic dipolar protic solvent and acts as a structure former. Hence, it 
removes water from the coordination sphere of metal ions, making them more 
reactive towards the ligands. As a result, the stability of the complexes is 
expected to increase. At the same time, it is a coordinating solvent and it com- 
petes with the ligands for coordinating the metals. This decreases stability of the 
complexes. Hence, stability of the complexes may increase or decrease. The vari- 
ation of overall stability constants with co-solvent content depends upon electro- 
static and non-electrostatic factors. Born's classical treatment holds good in 
accounting for the electrostatic contribution to the free energy change.” Accord- 
ing to this treatment, the energy of electrostatic interaction is related to dielectric 
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constant. Hence, the log B values should vary linearly as a function of reciprocal 
of the dielectric constant (1/D) of the medium. The linear variation observed in 
the present study (Figure 1) indicates that electrostatic forces are dominating the 
equilibrium process under the present experimental conditions. This is due to 
protophilic dipolar protic and structure forming nature of the solvent. The stabili- 
ties of the complexes follow the Irving-Williams order, i.e., Co(II) < Ni(II) << 
Cu(II). 


3.3. Effect of Systematic Errors on Stability of Ternary Complexes: 

In order to rely upon the best-fit chemical model for critical evaluation and appli- 
cation under varied experimental conditions with different accuracies of data 
acquisition, an investigation was undertaken by introducing pessimistic errors in 
the influential parameters like concentrations of alkali, mineral acid, ligand, 
metal, (Table 5). The order of the ingredients that influence the magnitudes of sta- 
bility constants due to incorporation of errors is alkali > acid > Asp > en > metal. 
Some species were even rejected when errors were introduced in the concentra- 
tions. The rejection of some species and increased standard deviations in the sta- 
bility constants on introduction of errors confirm the suitability of the experi- 
mental conditions (concentrations of ingredients) and choice of the best-fit mod- 
els. 


3.4. Stability of ternary complexes: 

The change in the stability of the ternary complexes as compared to their binary 
analogues was quantified” based on the disproportionation constant (log X) 
given by Equation | which corresponds to the equilibrium ML,+Mx, 
=—2MLX. Under the equilibrium conditions one can expect the formation of 
50% ternary complexes and 25% each of the binary complexes statistically and 
the value of log X shall be 0.6. A value greater than this accounts for the extra sta- 
bility of MLX. 
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Another approach” to quantify the extra stability of ternary complexes was 
based on the difference in stability (A log K) for the reactions ML with X and M 
(aq) with L and X, where L is the primary ligand (Asp) and X is the secondary 
ligand (en). It is compared with that calculated purely on statistical grounds as 
given in Equation 2. The equations for the calculation of A log K and log X are 
given in Table 6. 


AlogK =logK My a avd a cace (2) 


The electrostatic theory of binary complex formation and statistical arguments 
suggest the availability of additional coordination positions of the hydrated 
metal ion for the first ligand than for the second. Hence, the usual order of stabil- 
ity K“\,, > K™\, applies. This suggests that A log K should be negative, although 
several exceptions have been found.” The statistical values of A log K for 
tridentate L and bidentate X are, -6.67 to 7.06 in PG-water mixtures (Table 6). 
Negative values of A log K can be understood as the secondary ligand forms a 
more stable complex with hydrated metal ion than with L. Whenever the experi- 
mental values of A log K exceed the statistical values, it can be inferred that the 
ternary complex is formed as a result of interaction of ML with X or MX with L. 
The log X and A log K values could not be calculated for some systems due to the 
absence of relevant binary species. The log X values range from 2.12 to 18.66, all 
the values are found to be higher than those expected on statistical bases (0.6). 
These higher values account for the extra stability of the ternary complexes. The 
extra stability of ternary complexes makes them more amenable for metal trans- 
port. The less stable binary complexes make the metals bioavailable. 


3.5. Distribution diagrams: 

Based on the active forms of the ligands” (Lh, , LH,, LH, L” and XH,”, XH’ and 
X) and the existence of the binary complex species'’’'“(ML, ML,, ML;H,,ML,H, 
andML,H, and MX, MX, and MX,), the species are assumed to form through the 
following equilibria. 


ML,H,+ XH, — = MLXH+3H (1) 
M(Il) +LH+ XH, —= MLXH+2H" (2) 
M(I}) + LH + XH, MLX + 3H (3) 
ML + XH, —= MLX+2H’ (4) 
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MLXH —= MLX+H (5) 
ML.H,+ XH, — ML.X+4H (6) 
ML + LH + XH, ML,X + 3H’ (7) 
MLX + LH — MLX+H (8) 


The distribution of metal ions in various complex species (chemical speciation) 
with pH is represented as distribution diagrams. Some typical distribution dia- 
grams are given in Figure 2. The protonated ligands interact with metal ion 
(Equilibria 2, 3 and 4) to form MLXH and MLX. MLX is also formed by the dis- 
sociation of MLXH (Equilibrium 5) which is formed by the interaction of ML,H, 
with XH, (Equilibrium 1). Similarly ML,X is formed by the interaction of ML,H, 
with XH, (Equilibrium 6), ML and LH with XH, (Equilibrium 7) and MLX with 
LH (Equilibrium 8), because the concentrations of both MLX and LH are 
decreasing with increasing concentration of ML,X in the pH range 8.0-10.0. 


4. STRUCTURES OF COMPLEXES: 

The literature” suggests that Co(II), Ni(II) and Cu(II) complexes shall be octa- 
hedral. Amino nitrogen can associate with hydrogen ions in physiological pH 
ranges. Hence, there is often significant competition between hydrogen and 
metal ion for this donor site. This situation results in the simultaneous existence 
of a number of equilibria producing an array of protonated complexes, which are 
detected in the present study. Asps acts as bidentate or tridentate ligand depend- 
ing on the experimental conditions and en acts as bidentate ligand. Thus based on 
the above equilibria the speculative structures of the complexes are presented in 
Figure 3. 


5. CONCLUSIONS: 

The following conclusions have been drawn from the modeling studies of the 
speciation of ternary complexes of Co(II), Ni(ID) and Cu(II) with Asp and en in 
PG-water mixtures. 


1. The ternary metal complex species detected are MLX, ML,X and 
MLXH for Co(II), Ni(I]) and Cu(II), where L=Asp and X= en. 


2. The values of A log K and log X indicate that the ternary species have 
extra stability compared to their binary species. 


3. The linear increase in the stabilities of ternary complexes with decreas- 
ing dielectric constants is due to the dominance of electrostatic forces 
over non-electrostatic forces. 


4. The Alog K method is based on statistical grounds and the log X method 
takes dispropotionation of the binary complexes into consideration. A 
log K should be negative if the ternary complex has extra stability. Any 
value for log X greater than 0.6 accounts for extra stability. 


5. The order of the ingredients that influence the magnitudes of stability 
constants due to incorporation of errors is alkali > acid > Asp > en> 
metal. 


6. The study also gives an insight into the metal availability/metal trans- 
port in biofluids. The ternary complexes are more amenable for “metal 
transport” because of their extra stability. 
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PG 


Temperature = 303.0+0.1 K; Mineral acid = 1.0 mmol; p = 0.16 mol L”. 


Tm0 


Table 1: Total initial concentrations of reactants (in mmol) for mixed ligand titrations in PG-water mixture. [NaOH] = 0.4 mol L”; V, = 50 mL; 


. - TM0:TLO:TX0 
% VIV Co(II) Ni(II) CudD 
0.2500 0.2499 1D S95 
0.0 0.100 0.099 0.102 0.2500 0.5000 1:2.5:5.0 
0.4999 0.2499 1:5.0:2.5 
0.2498 0.2499 IDSA 
10.0 0.100 0.099 0.102 0.2498 0.5000 1:2.5:5.0 
0.4998 0.2499 1:5.0:2.5 
0.2492 0.2500 1255 
20.0 0.100 0.099 0.102 0.2492 0.5000 1:2.5:5.0 
0.5001 0.2500 e525 
0.2510 0.2502 DSS 
30.0 0.100 0.099 0.102 0.2510 0.5000 1253.0) 
0.4998 0.2502 1253,.022.5 
0.2550 0.2588 1D 50.5 
40.0 0.100 0.099 0.102 0.2550 0.5000 1:2.5:5.0 
0.4999 0.2515 1:5.0:2.5 
0.2500 0.2505 1D 52S 
50.0 0.100 0.099 0.102 0.2500 0.5000 1255.0) 
0.5001 0.2505 1:5.0:2.5 
0.2500 0.2503 DBD. 
60.0 0.100 0.099 0.102 0.2500 0.5000 1D SS.0) 
0.5001 0.2503 1:5.0:2.5 


Model No. 


log B(SD) 
ML,X 


303.0+0.1 K, Ionic strength = 0.16 mol L”. 


Skew-ness 


Kurto-sis 


Table 2: Exhaustive modeling study performed on Asp-Ni(II)-en ternary complexes in aqueous medium, pH range =4.0-10.0, NP=92; Temperature= 


R-Factor 


1 11.96(12) : z 3.84 3.46 0.81 4.82 0.0008 
5 S 22.18(15) z 3.44 | 5.32 0.20 10.65 0.0012 
3 ; : 27.58(7) 20 ae 76.96 0.0015 
4 12.91(12) 21.30(20) = 2.04 | 1.15 -1.06 2.43 0.0109 
5 13.56(8) : 24.33(6) 1.80 le 0.97 4.26 0.0006 
6 = 20.73(16) 22.12(8) 1.69 | 3.42 1.29 28.91 0.0075 
7 15.01(32) 19.87(17) 21.15(10) 0.84 naan 23.75 773. 00177 


Table 3: Best fit chemical models of ternary complexes of Co(II), Ni(II) and Cu(II) with Asp and en in PG-water mixtures. 


1111, 


Log Byus(SD) 
1110 


1210 


NP 


Co(II) (pH=5.0-10.0) 


U,,.X10* 


corr 


2 


x 


Skewness 


Kurtosis 


R-factor 


00.0 19.33(12) 1.888) 15.5123) 23 0.85 12.35 -0.37 3.65 0039 
10.0 19.46(12) | 12.15(12) | :16.16(12) 37 1.52 116.4 2.33 7.97 0109 
20.0 19.67(22) | 12.64(28) | ——-16.7(20) 35 2.68 69.28 3.07 12.75 0289 
30.0 19.87(12) 12.8312) «16.86 (12) 48 1.90 74.00 3.46 17.27 0526 
40.0 19.99(12) / -12.86(12) | ——«16.93(12) 32 0.16 34.88 1.10 3.71 0432 
50.0 20.18(12) Sea aie a ee [4.14 12.94 2.01 8.80 0064 
60.0 20.78(12) | 13.6412) | —«:18.37(12) 23 1.52 12.36 1.75 7.56 0132 
Ni(ID (pH=4.0-10.0) 
00.0 21.15(57) 15.0182) 19.8733) 9% | 7.84 | 23.75 1.60 7.23 0178 
10.0 21.95(70) 15.3949) «19.97(42) 34 9.38 22.09 2.12 8.57 0163 
20.0 22.28(24) [| 15.5922) | 200139) | 36 | 6.18 14.55 -0.15 3.01 0134 
30.0 22.33(25) | 15.72(28) | —-20.24(36) 35 3.14 32.33 2.47 11.20 0094 
40.0 22.6018) | 16.0205) | 204009) | 25 «| 3.28 4.32 -0.22 3.09 0093 
50.0 22.61(11) | 16.33(12) | —_20.71(15) 26 0.60 23.08 2.79 15.08 0041 
60.0 23.27(32) | 17.1163) | 22.6054) | 27 9.54 12.59 2.61 10.07 0178 
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Cu(ID (pH=2.0-8.0) 
00.0 30.31(24) 25.56(45) |  32.65(27) 30 2.66 5.87 -0.19 3.36 0096 
10.0 30.09(10) | 25.61(25) | 33.12(22) SH/ 32D) 4.63 -0.36 3.07 0085 
20.0 30.83(30) | Pee) | 33.43(33) 33 2.83 16.11 -0.42 3.38 0117 
30.0 31.00(42) | 26.28(49) | 33.58(47) ae) 4.58 48.07 -0.38 3.58 0199 
40.0 31.22(14) | 26.49(40) 34.88(23) 36 1.07 32.44 2.01 10.93 .0083 
50.0 31.38(13) | 27.66(13) 35.01(13) 37/ 0.20 Tiel -0.06 4.84 0016 
60.0 32.24(4) | 28.18(8) | 35.24(7) | 39 0.95 10.15 -0.46 3.09 .0034 
Table 4: Comparison of ternary stability constants reported in literature with the present study. Instrumental method is pH metry. 


log Buus 


—--— Sener aati itaaaenealeaaaaaaate Ionic strengthh (M) Reference 
1110 1210 1111 1120 1121 
Co(II)- (L-Asp)-(en) iL S7/ 14.17 | - 14.83 24.29 S1ES2 0.16 29 
Co(II)- (L-Asp)-(en) 11.88 IS.s511 I@).3%3 - - - 0.16 Present study 
Ni(II)- (L-Asp)-(en) 13.89 16.92 - 18.35 26.62 | 32.46 0.16 29 
Ni(II)- (L-Asp)-(en) 15.01 IO | QL IS - - - 0.16 Present study 
Cu(ID- (L-Asp)-(en) 17.81 20.98 - AILEY) 30.40 37.11 0.16 29 
Cu(ID- (L-Asp)-(en) 25.56 32.65 | 30.31 - - - 0.16 Present study 


Table 5: Effect of errors in influential parameters on the stability constants of Asp-Co(I])-en ternary complexes in 50% v/v of PG-water mixture. 


i xn(SD 
Ingredient 08 Baw(SD) 
1110 
0 20.18(12) 13.34(12) 17.10(12) 
-5 Rejected 13.71(29) Rejected 
ek -2 16.23(27) 15.70(38) 19.11(48) 
ali 
+2 19.92(41) 15.00(35) Rejected 
35) Rejected Rejected 19.15(29) 
-5 19.31(16) 12.82(29) 18.33(29) 
ee -2 18.90(23) Rejected 18.96(20) 
ci 
42) Rejected 14.71(35) | 14.82(99) 
+5 18.35(43) Rejected Rejected 
-5 18.73(35) 14.52(67) | 16.30(30) 
rn -2 18.25(26) 14.17(46) 17.85(27) 
423) 18.20(25) 14.02(42) 16.15(36) 
+5 17.97(28) aL 13.98(38) 15.83(79) 
-5 17.56(30) 13.86(17) 17.74(42) 
-2 18.12(26) 13.97(35) | 17.23(32) 
en 7 
+2 18.59(41) 13.98(78) | 17.82(43) 
+55) 19.98(38) 13.69(58) | 17.29(40) 
-5 19.55(30) 14.20(53) | 18.99(28) 
Mea -2 19.43(29) 14.04(48) | 18.02(31) 
eta 
+2 19.27(54) 14.06(57) 18.14(62) 
+5 19.14(24) 14.64(37) | 18.59(39) 


Table 6: A log K and log X values of ternary complexes of Co(II), Ni(II) and Asp-Cu(II)-en in PG-water mixtures. 


meh A log K log X 
Co(II) 
1110 1210 1110 | 1111 i 
0.0 -0.47 5.77 2.48 2.71 
10.0 -0.39 -5.75 2.60 BOT 
20.0 -0.37 -6.38 2.20 2.73 
30.0 -0.29 -5.75 3.05 3.15 
40.0 -0.51 -6.67 2.12 3.69 
50.0 -0.03 -6.50 2.86 5.05 
60.0 -0.19 -6.22 2.69 4.67 
Nid) ii 
0.0 0.13 -5.19 3.68 4.69 
10.0 -0.58 5,56 3.27 4.40 
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20.0 -0.33 5.12 3.78 4.16 
30.0 | -0.22 | -5.49 | 3.62 | 5.65 
40.0 -0.41 | -6.24 | 3.20 | 5.69 
50.0 -0.34 Paes ay) 6.12 
60.0 | 0.51 -6.08 4.87 7.64 
Cu(II) 

0.0 ae BA 15.23 15.77 
10.0 6.28 -2.49 15.61 15.78 
20.0 | 5.21 -4.93 13.14 5 

30.0 4.96 -4.50 14. 48 16.14 
40.0 5.38 DE 15.34 16.46 
50.0 | 6.21 BD | 17.09 | 15.51 
60.0 7.06 -2.46 18.66 17.55 


Equations for the calculation of A log K and log X 


A log Ku =log Binto -log Buroo -log Broo 
A log Kuno =log Bisto -log Bisoo -log Bioxe 
A log Kun =log Bun -log Baroo -log Brow 
log Xin =2log Brno -log Bisoo -log Broo 
log Xin =2log Biaw -log Braoo -log Broan 
log Xun =2log Bun -log Biooe -log Broan 
m B 
24 
20 
22 
18 
es O° = »| gee 8 - 
16 t=?) 
e 2 
_ 18 
14 
16 
12 
0.012 0.016 0.020 0.012 0.016 0.020 
4/D 1/D 
36 
Cc 
33 
= 
wD 30 
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27 
0.012 0.016 0.020 
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Figure 1: Variation of stability constants of Asp-metal-en ternary complexes in PG-water mixtures. (A) Co(ID, (B) Ni(II and (C) Cu(ID (0) log Burx. 
(0) log Byuox and (A) log Bonixe. 
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Figure 2: Species distribution diagrams of ternary complexes of Asp and en in 40% v/v PG-water mixtures. (A) Co(ID, (B) Ni(ID and (C) Cu(II). The 
amounts of ingrediants are M(II): 1.0 mmol, Asp: 5.0 mmol and en: 2.5 mmol. 
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Figure 3: Speculative structures of Asp-M(II)-en complexes, where S is either solvent or water molecule. 
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